Zero-point energies and the multiplicative anomaly. 
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For the case of a relativistic scalar field at finite temperature with a chemical potential, we 
calculate an exact expression for the one-loop effective action using the full fourth order determinant 
and ^-function regularisation. We find that it agrees with the exact expression for the factored 
operator and thus there appears to be no mulitplicative anomaly. The appearance of the anomaly 
for the fourth order operator in the high temperature limit is explained and we show that the 
multiplicative anomaly can be calculated as the difference between two (^-regularised zero-point 
energies. This difference is a result of using a charge operator in the Hamiltonian which has not 
been normal ordered. 
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I. INTRODUCTION 



Renormalisation and regularisation techniques form a vital part of the physicist's arsenal when performing cal- 
culations in quantum field theory. The technique of ^-function regularisation jl]-[| is a well established method for 
obtaining finite results in quantum field theoretical calculations and has shown itself to be a very elegant and powerful 
method, particularly in curved spacetimes and spacetimes with a non-trivial topology. 

Calculations involving Feynman path integrals typically involve the determinant of a differential operator. This 
determinant is an infinite product and has to be regularised in some way. It has been shown recently that results 
04 • obtained from ^-function regularised determinants can be ambiguous - the source of this ambiguity is known as a 
multiplicative anomaly. The multiplicative anomaly is essentially the difference between different zeta regularised 
i factorisations of a determinant. For example, in calculating a determinant one may wish to factorise it in order to 
make calculations easier. Normally one would write det(AB) = det(A) det(-B), but in the case of infinite matrices 
this relation is not always correct after regularisation. The multiplicative anomaly in a D dimensional spacetime is 
i-H ■ defined as, 



Qh! a D = lndet(AB) - lndet(A) - lndet(B) (1.1) 

CD 

■ The relevance of the multiplicative anomaly for physics was first brought to light by Elizalde, Vanzo and Zerbini 

and they showed its connection with the Wodzicki residue. The high temperature limit of the one loop effective 

■ action for a charged scalar field with chemical potential was considered by Elizalde, Filippi, Vanzo and Zerbini ||, 
and the resulting anomaly was found to depend on the chemical potential. It appeared therefore that an extra term, 
previously overlooked, might be present in the effective action, a term which could not be removed by renormalisation. 
(This idea received criticism from Evans and Dowker p). Elizalde, Filippi, Vanzo and Zerbini responded to these 
criticisms in |||l0[.) 

It seemed that there were many different expressions for the effective action, one for each way in which the deter- 
minant can be factorised. So, given these varying expressions, each one differing from another by a corresponding 
multiplicative anomaly, how can we know which one (if any) is correct? The present authors concluded in [|ll| that 
the ambiguity associated in choosing different factorisations could only be resolved by making comparisons with cal- 
culations performed using canonical methods. In this paper we shall consider two factorisations; the full fourth order 
operator: 



1 ■ - *4 ft n i ™2 ,,2\ 2 a ,,2 ^ 



r A = -lndeU 4 -D+m J V "V^To (1-2) 



2 ^ ■ f , dt 2 

and the 'standard' factorisation of two second order operators (eg. as used in |12fl): 
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(1.3) 



In [jrij an exact expression for the effective action could only be calculated for the B factorisation, and this agreed 
with the standard, well known thermodynamical expression - a sum over the zero-point energies and the thermal 
contributions of Bose-Einstein sums for particles and anti-particles. The high temperature limits of the effective 
action for both cases were calculated and it was found that T b agreed with the result obtained by Haber and Weldon 
p3[ who did not use path integrals or C- function regularisation. Ta differed from Tb by an amount exactly equal to the 
multiplicative anomaly calculated in H . There does not seem to be an a priori way of determining which factorisation 
will yield the correct physics; an 'objective' comparison with canonical methods needs to be performed. Clearly this 
is a problem if calculations using ^-function regularisation need to be made on a system where the canonical answer 
is not known. 

Since the high temperature expansions of and Tb differ, it might be thought that there would be some discrepancy 
between their exact expressions also. In Sec. || we shall calculate the exact effective action Ta, and show that it 
actually gives the correct result, in complete agreement with T^. Is the multiplicative anomaly therefore just an 
artefact of the high temperature expansion? In Sec. Ill we postulate that the multiplicative anomaly arises from the 
zeta-regularised zero-point energies and is not a thermal phenomenon. It arises in factorisations where the charge 
operator Q has not been normal ordered. The multiplicative anomaly is calculated explicitly as the difference between 
^-regularised zero-point energies with and without a chemical potential. We also consider the interacting case. Both 
calculations give multiplicative anomalies which agree with those calculated in 0. In Sec. [V we shall draw some 
conclusions. 



II. THE EXACT EFFECTIVE ACTION F A 

We shall use the notation and conventions of now setting e = 1. We are working with a relativistic, non- 
interacting, charged scalar field with a chemical potential. The action is 

S[(j>] = J dt J da x - i^4>2) 2 + ^(4>2+ i^4>if 

+ ^|V0i| 2 + i|V0 2 | 2 + im 2 (^ + ^)^ . (2.1) 

We expand about a constant background field with <f>\ — (f> and <p2 = 0. T a is defined formally as 

1 



T A = -IndeU 2 (l 2 S^[ 



= -lndeU 4 f (-n + m 2 -/, 2 ) -4 M 2 ^). (2.2) 



Using ^-function regularisation || we can define: 

1 ... 1 



2 CA(0) + ^a(0)ln/ 4 (2.3) 



with 

CO _ 

Us) = E E [K - ^ + Elf + 4^ 2 ] " (2.4) 



n J — — QO 



The ujj are the Matsubara frequencies for scalars, cuj — 2irj / f3 and = & n + m 2 . a n are the eigenvalues of —V 2 on 
the spatial part of the manifold. We shall apply the Abel-Plana summation formula, 

°° /*+oo roo-\-ie 

£/(i)=/ f(x)dx+ dz(e-^ z -iy 1 [f(z) + f^z)} (2.5) 

_ ~~ J —oo J— oo+ie 



j = -oo 



to evaluate Ca(s). Let us label the two integrals arising on the right hand side of ( |2.4| ) Q(s) and R(s) respectively 
after using (|2.5|) to perform the sum over j. Then, 



2 



= / [ax 4 + frx 



+ 



dx 



(2.6) 



where a = (2tt/0) 4 , b = 2(2n / (3) 2 (E 2 , + fi 2 ) and c = (J5{j — ^ 2 ) 2 . Making a simple substitution and using the identity 

1 



a 



gives us 



(fox-* e -(a* 2 +bx)t _ 



Evaluating the integral over x first results in 

1 



Q00 = £ 



d£ 



" V 8a 



f-l if! I ^ 



2r(s) V a J 

where Ki is a modified Bessel function. Performing the integral over t |Q we get, 

1 3 1 



(2.7) 



(2.8) 



(2.9) 



n( ,_sr 1 /I ^ /fr 2 yr(.s + i)rQs-i 

WJ Z^ 2 r( S )Vac s + 1 /4 ^4a 



r( s + i) 

4- V 2 ^2^; _ M 2 )25+ i/ 2 r(s)r( i) y o ^ 



2-F1 s + -, -:s + -; 1 
" 4 4 2 



-1 

4ac / 



(2.10) 



after substituting in the values of a, b and c. This is analytic at s — 0, and is equal to zero, Q{0) = 0. It is easy to 
see that, 



n * 2 1>J (i? 2 -, 2 ) 1 / 2 r ( i) ^ixl) I 1 + 

n 

= -2PY, E n 



{El - M 2 



1 

— arctan 
2 



{El-H 2 ? 



arctan 



2^E n 

{ E n — M 2 ) 



(2.11) 



This last result involves the sum of zero-point energies. Turning now to the contour integral in (2.5), we can write 



-2ttiz 



— oc+ie 



2-kz 



- + El 



4 M 2 



2?rz 

T 



dz 



(2.12) 



since /(z) is even. There are poles at all integers on the real axis, and branch points where the expression in square 
brackets in the above equation is equal to zero, namely at z — ±i(/3/27r)(i?„ + /i),±i(/3/27r)(_E„ — /i). By taking branch 
cuts between the poles in each half of the complex plane, and closing the contour in the upper half plane, it can be 
shown that 



R(s) = -2s£ln [(l - e-^"-")) (l - e~^ + ^ 



Hence R(0) = and 



#(0) = -2£m[ 



1 - e 



-/3(B„-At) 



)( 



1 - e 



(2.13) 



(2.14) 
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Therefore £a(0) = Q(0) + i?(0) = and so 



Ta = -^(0) 



^ + In [(l - e -^ E ^) (l - e -^ E ^)] } (2.15) 



which is the same result as was obtained for Tg in |ll|. A simpler, but perhaps less elegant way of evaluating (2.4) 
is presented in the appendix. 

So we have a paradox - the exact expressions for Ta and Tb agree, while their high temperature expansions differ 
by the multiplicative anomaly. The resolution of this paradox lies in the fact that the sums over the energy levels 
(integrals over k when a n — k 2 ) have not yet been performed in the exact expressions we are considering. In the high 
temperature expansions, the ^-functions were expanded in powers of /i, the chemical potential, and the sums over the 
energy levels and Matsubara frequencies were then performed. As we shall see in the next section, the multiplicative 
anomaly arises from the chemical potential being present in the zero-point energy contributions. This is due to a lack 
of normal ordering in the charge operator. 

III. NORMAL ORDERING AND THE MULTIPLICATIVE ANOMALY 

In this section we shall show that the multiplicative anomaly stems from the zero-point energy contribution to the 
effective action. 



A. Non-interacting model 

There are two ways to write down the zero-point energies for the system described by ( [2. if ); with or without a 
chemical potential. One can write {(3/2) ^2 n {E n ±/i) for particles and anti-particles (— /x and +/i respectively) or one 
can simply write (3^2 n E n (which was derived in our exact expressions for Ta and Tb)- The multiplicative anomaly 
is the difference between these (^-regularised zero-point energies. Let us define 

n 
n 

Formally of course, there is no difference between the zero-point energies in the two cases: i+ + i_ — J = 0. But if 
we regularise first then the anomaly appears. Define 

7 ±( s ) = f £(£n±M)-* (3.3) 

n 

J(s)=pJ2 E n S - (3-4) 

n 

Then, we claim that the multiplicative anomaly is 

04 = /+(-l) + J_(-l)-J(-l). (3.5) 
Let us calculate I±(s) with E 2 = k 2 + m 2 : 

I±(s) = f ^ J q °° dk k 2 [(fc 2 + m 2 Y± M ] ~ S (3.6) 

We can binomially expand the square bracket in powers of \i up to 0(/x 4 ). We do not need to consider higher order 
terms for reasons which will become apparent in due course. After performing the integrals we have, 
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Ts(s + i )(s + 2)^^^!^ (m 2 )- f + + l)(.s + 2)( S + 3)fJ^| (m 2 )"^ j . (3.7) 

The terms with odd powers of // cancel when we write down an expression for i+(s) + i_ (s). In the fi 2 and /x 4 terms, 
the T- function in the numerator is divergent at s = — 1, but can be analytically continued to cancel away the factor 
of (s + 1) multiplying each term. Thus, 

™ + '-<•> " ^ { V M 1 " 5 + ^.i 2 ^ 

+s(s + 2)( S + 3)^^|±^ H -^| . (3.8) 

All even, higher order terms in /i have analytic T-functions in the numerator, and so the (s + 1) ensures they are all 
zero at s = —1. This is why we were able to stop expanding at f ourth order in /i. 
Turning now to J(s), we see it is exactly the first term of (3.S): 

J( s ) = l^ ^ dkk 2 (fc 2 +m 2 )" f 

(27Tj Jo 

V^r(l-I) 2 j_ f 

" 8vr2 r(f) lW j (d ' Jj 



Using (3.5), it is now a straightforward matter to show that the multiplicative anomaly is 

«4 = ^M^m- y j (3.10) 

in agreement with [||. 

This calculation sheds some light on why the high temperature expansions of Ta and Tb differ, and the exact 
expressions agree. In the high temperature situation, the integrations over k were carried out after the expansions, 
and for a reason which is not clear to us at the present, the chemical potentials in the zero-point energies of the 
A-factorisation were not able to cancel. So the zero-point energies were of the form {/3/2)^2 n (E n ± /u). In the 
B-factorisation, the energy levels were just (3 J2 n E n and so there was no anomaly. In the exact expressions for Ta 
and Tb, the integrals have not even been performed, and so the and —fi simply disappear, leaving no trace of a 
discrepancy. 

Although simply by looking at the A factorisation (2.2), (2.4), we cannot say whether or not it will produce a 
multiplicative anomaly in the high temperature expansion, given that we know it does produce an anomaly, we can 
say something about why it does. The canonical energy levels for our system are derived from the Hamiltonian 
operator H, 



H = Y, En(ala n + ±+ 6+6„ + ~) (3.11) 



where ajj, a n (6^, b n ) are the creation and annihilation operators for particles (anti-particles). For a system of 
charged fields with a chemical potential, the full Hamiltonian (which is the argument of the exponential in the 
partition function) is 

H = H-^:Q: (3.12) 
where : Q : is the normal ordered charge operator 

: Q := X>K - bib n ) . (3.13) 

n 

Note that we have to normal order by hand. There is no good mathematical reason why we normal order, we just like 
to have an uncharged vacuum: 
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(0| :Q: |0> =0 



(3.14) 



It now becomes clear why we have two different expressions for the energy levels ( |3.l|) , (p^). They correspond to 
the eigenvalues of H and H respectively in the case when the charge operator Q is not normal ordered. So to avoid 
having an anomaly we need to ensure that both H and H have the same eigenvalues - we need to normal order Q. 
This implies that the A factorisation, which gives rise to an anomaly, is not normal ordered. This is a symptom of 
using Fcynman path integrals, indeed Bernard was aware of this in 1974 - the last sentence in section II of his seminal 
paper |l5[ reads, '. . .the functional-integral formalism never does normal ordering for us.' 



B. The interacting case 

The multiplicative anomaly can also be calculated in the interacting case, and can again be seen to be the difference 
of ^-regularised zero-point energies. 
Let us define 



E 



El 



\<t> 2 



A0 2 \ AW 



3G 



2 \ 2 



6 



2 \ 2 



(3.15) 



(3.16) 



in an analogous way to I± and J in the non-interacting case. (See for example fll6fl for a full derivation of the energy 
levels). Then we have 



El 



n*) = ?EK 



\<t> 2 



M ± V E, 



3(» 







EK 



A^ 
6 



2J,4\ 2 



(3.17) 
(3.18) 



To evaluate X±(s), we expand the square root inside the square bracket in powers of A up to 0(A 2 ): 



F 2 , V , a , _ „ , mA0 2 A 2 4 M A 2 ^ 4 



3£„ 144 M £„ 36££ 



(3.19) 



We note that we can switch from X+(s) to X-(s) by letting /j, — > — Therefore we shall work with -X+(s) for 
simplicity, and let /1 — > — /i to write down X_(s) at the end of the calculation. X + (s) can be written in the form, 



X+(s) = | 2 [ao + 01 A + a 2 A 2 ] 



21-2 



a 2 [l + a 1 aiA + a ^A 2 ] 



I EK 



1 



1 



2 „2 \2 



— sa 2 aiA sa 2 a2A H — s(s + 2)a 2 a^A 



(3.20) 



where 



a = (S„ + /1) 

°i = y (i + z^n 1 ) 

° 2 - 36 [nr - 



(3.21) 



It is then a straightforward (if a little tedious) matter to expand each term of (3.20) and integrate over k, as was done 
in the non-interacting case. 
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After the dust settles, we find 



V^pv (r(§-§) 



,2\ 2 _ 2 



^ + x -^-^i^f(i) — — (s -i)r V2 



2s/i 



2 r(f + |) 



T) ' ' 



2\ 2 _ 2 



+s(s + 2)(s + 3) 



M 4 r(§ + |) 2 ,-|-U ^/3va</> 2 ( r(§ - ±) 



12T(§ 



48tt 2 



r(§ + i) 



fy(™ 2 ) 



2\ 2 _ 2 



(to 2 )" 



5v^/3VA 2 (/- 4 s(s + 2) T(i 



2304tt 2 



r(f + 2) 



(to 2 ) 



2\ - 2 _ 2 



By expanding ( 3.18| ) to order A 2 , it can easily be shown that 

I) 

Cto 2 V 

6 r(f + i) [ > 



^) = ^(M M l-f_ s ^ r (f-i) 



stt 2 l rm 



r(§ + 2) 



(3.22) 



(3.23) 



and hence, 



a 4 = X+(-l)+X_(-l)-y(-l) 
/?F 2 / 2 A0 2 /i 2 

= 8^ m + — "T 



(3.24) 



This agrees with the result in ^] except for a term in A 2 4 . This is of no consequence however, since any term 
proportional to the background field 4> (but not the chemical potential) may be added to the effective action without 
changing the physics of the system. All such terms can be harmlessly absorbed by renormalisation. 

This section has shown that the multiplicative anomaly has its roots in the manipulation of infinite sums - the 
non-interacting case in particular demonstrates how the anomaly can appear in relatively simple situations. Elizalde 
showed the existence of the multiplicative anomaly in p ossibly the simplest of all cases - infinite, diagonal matrices 
with real numbers |]l7f . The first worked example in |l7f | is striking in its similarity to the calculation presented above 
in the non-interacting model. 



IV. CONCLUSIONS 



We have demonstrated that the multiplicative anomaly originates in the zero-point energies of fields and is a 
consequence of shifting the energies by a constant amount. When regularisation is performed, these shifts (+/x and 
— n for example) are unable to cancel and result in a multiplicative anomaly. It seems therefore that in order to 
perform anomaly-free calculations one must resist the temptation to integrate over the momentum until the very end, 
after say, a high temperature expansion has been written down. 

It should be borne in mind that the functional integral approach to quantum field theory is not as complete as the 
canonical one. As was mentioned in Sec. Ill, path integrals completely neglect normal ordering. Coleman discusses 
the merits of path integrals at length in his Erice lectures and echoes the comments of Bernard; the functional 
integral approach does not norm al or der. 



From inspection of equations (3.1) and (3.2) it is tempting to conclude that C~ function regularisation may be to 
blame for the multiplicative anomaly (as opposed to functional integration). Certainly it does not see m lik ely th at a n 
anomaly would survive if say, dimensional regularisation were used to calculate the difference between (3.1) and ( [3.2] ). 
But it should be remembered that these equations were written down almost naively, to demonstrate the source of the 
anomaly; as was mentioned above, one should wait until the last possible moment before performing the integration 
over k, after everythi ng t hat ca n ca ncel has done so. Nevertheless this does not seem very satisfactory. Why is there 
a difference between (3.3) and (3.4)? The mathematical properties of the ^-function are rigorous and well defined - 
to negatively criticise the whole subject of ^-function regularisation is a step not to be taken lightly. The paper by 
Elizalde [jlTj provides some very interesting mathematical examples of multiplicative anomalies derived from infinite 
matrices, sometimes using nothing more than Riemann's ^-function. 

These mathematical peculiarities aside, it has been clearly demonstrated in this paper that the problem associated 
with the multiplicative anomaly can be removed by considering a Hamiltonian with a normal ordered charge operator. 
Normal ordering is in some ways an artificial procedure that physicists perform to make the theory more physical - 
it is not prescribed by the theory and there is no mathematical reason why it is done. Using canonical techniques it 
is easy to see how to normal order, unfortunately it is not so obvious in the functional integral approach - equation 
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(1.2) is not normal ordered, but ([0]) is. An important step in understanding the multiplicative anomaly would be 
finding an a priori method of knowing which factorisations lead to anomalies and which do not, without having to 
compare with canonical calculations. 

The problem of the multiplicative anomaly appears to be quite deeply rooted in the mathematics of infinite, 
divergent series. 
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APPENDIX: AN ALTERNATIVE METHOD OF CALCULATING V 



(i) 

A 



The method presented here for calculating (3.5) is simpler than that given in the main text, but has the disadvantage 
that the function can only b e eva luated at s = 0. The method in the main text can be used for an arbitrary value 
of s (although the i nteg ral in ( 2.10 ) may have to evaluated numerically). 

We can re-write (2.4) as, 



0i(s;a,6) 



E E (f+-r s (f+b 2 y 



n j — — oo 

where a = (3/2ir(E n + /j,) and b = (3/2n(E n — fi). Differentiating with respect to (a 2 ) and (b 2 ): 

d 

— t^-Ca(s; a, b) = -sf(s; a, b) 
o(a z ) 

d 

q^(a{s; a, b) = -sg{s; a, b) 

where 



(Al) 



(A2) 



Clearly, 



f{s;a,b) 
g(s;a,b) 



2 i) EE 



(3 



n J — — OC 
-As oo 



2 EE 



(/+« 2 ) 



(f+b 2 )- 8 ' 1 



n j — — oc 



d 



-U(0;a,b) = 



and so we can conclude 



d{a 2 Y 
d 

qK 2jUwa,b) = 



U(0;a,b)=C 



(A3) 



(A4) 



(A5) 



where C is a constant independent of a and b. We can set a = b in order to evaluate the left hand side of (A5) and 
determine the constant C. Using the Plana formula to calculate Ca{s] a, a) is straightforward: 



°° poo noo+ie 

E (i 2 +« 2 )" 2S =/ dj(f + aY 2S + 2 (e-^-irV + ^T 2 ^ 

- J — oo J — oo-t-ie 



(A6) 



and so 



C,A{s;a,a) 



(A7) 
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(See the appendix of |Tl| for a more thorough evaluation of a similar sum.) Consequently £a(0;<i, a) = C = 0. Next 
we need to evaluate £' A (0;a,b), 



d 



d(a 2 
d 

d(b 2 



■C' A (0;a,b) = -f(0;a,b) 
-C' A (0;a,b) = -g(0;a,b) 



The functions /(s; a, b) and g(s; a, b) can easily be calculated at s — 0: 
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d{a 2 ) 

d 
9(62) 



Ca(0j °j — — — coth(7ra) 

n 

C;(0;a,6) = -^Jcoth(^6). 



Integrating, 



C A (0; a, b) = {-2 In [sinh(Tra)] - 2 In [sinh(Trfe)]} + K. 



Again there is a constant IC, independent of a and b. We can evaluate it in the same way as before, 

K, = C' A (0;a,a) + 4 ^ In [sinh(Tra)] . 



Using (A7) we see that K. = —4 In 2. Writing the sinhs in ( |A10| ) as exponentials we arrive at 
C A (0;a,b) =^{-27ra-27rfe-21n(l-e- 27rQ ) - 2 In (l - e" 27rb ) } . 



Substituting in the values of a and b, 

Ca(0) = { + M) - P( E n ~ M) ~ 21n (l - 6-^"+^) - 2 In (l - e 

n 

and so finally we have, 



)} 



r 



(i) 



;&(0) 



Y {pE n + In (l - e-n E » + ri) + In (l - e -/*(*.-*0)} 



(A8) 



(A9) 



(A10) 



(All) 



(A12) 



(A13) 



(A14) 



It is interesting to see that the zero-point energies in (A13) are written initially with positive and negative chemical 
potentials. 
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